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We consider the signatures of a population of primordial black holes (PBHs) in future observations 
of 21cm radiation from neutral hydrogen at high redshift. We focus on PBHs in the mass range 
5 X lO^^kg < MpBH ^ lO^^kg, which primarily influence the intergalactic medium (IGM) by heating 
from direct Hawking radiation. Our computation takes into account the black hole graybody factors 
and the detailed energy dependence of photon and absorption by the IGM. We find that for black 
holes with initial masses between 5 x lO^^kg < A/pbh ;$ lO^^kg, the signal mimics that of a decaying 
dark matter species. For black holes in the range 5 x 10^"kg < A/pbh J; 5 x lO^^kg, the late stages 
of evaporation produce a characteristic feature in the 21cm brightness temperature that provides a 
unique signature of the black hole population. If no signal is observed, then 21cm observations will 
provide significantly better constraints on PBHs in the mass range 5 x lO^^kg < Mpbh ^ lO^'^kg 
than are currently available from the diffuse 7-ray background. 



I. INTRODUCTION 

The 21cm hyperfine spin-flip transition of neutral hy- 
drogen (HI) may allow observers to probe the cosmic 
"Dark Ages" comprising the epoch between the last scat- 
tering of the cosmic microwave background (CMB) at 
z 1100, and the appearance of luminous sources at 
z - 30 ^, 

0,3- Observations at wavelengths 21(1 -I- z) 
cm could be used to slice the universe as a function of red- 
shift z, and so produce a three-dimensional map of the HI 
distribution when combined with angular measurements. 
At moderate redshifts, this could allow reionization of the 
universe to be studied in detail. Overall, the essentially 
three-dimensional nature of 21cm data gives it the po- 
tential to become one of the the richest data sets for cos- 
mology [4]. It will not only enable the standard ACDM 
cosmological model to be better tested and understood, 
but will provide a new means to constrain - or detect - 
more exotic possibilities. One such possibility, the pres- 
ence of a population of primordial black holes (PBHs), is 
ideally suited to study with the 21cm background, as we 
describe in this work. 

The 21cm signal is sensitive to exotic physics through 
its dependence on the thermal history of the intergalac- 
tic medium (IGM). The IGM is visible in 21cm when 
the spin temperature Tg of the neutral hydrogen gas dif- 
fers from the CMB temperature Tcmb- The spin tem- 
perature Ts is itself determined by the competition be- 
tween its coupling to Tcmb through interactions with 
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GMB photons, its coupling to the gas kinetic temper- 
ature Tk through atomic collisions, and (especially at 
lower redshifts) interactions with Ly-a photons produced 
by luminous sources. Any process that heats the IGM 
will influence Tk and therefore affect T5. For redshifts 
z - 30-300, one expects - lO-lO^X - lO^^-io-i 
eV, and so little heating per hydrogen atom is necessary 
to significantly change the thermal history of the IGM 
and influence the 21 cm signal. By contrast, the CMB 
is mainly sensitive to the Dark Age IGM thermal his- 
tory through the effects of IGM heating on the ionized 
fraction, which affects the optical depth to the last scat- 
tering surface tlss- The CMB is therefore less sensitive 
to IGM heating, since the energy required to ionize HI is 
much larger than the typical Tk over the redshift range 
of interest. Furthermore, changes in the CMB TT power 
spectrum due to tlss a-re degenerate with changes due 
to the scalar spectral index (though measurements of 
the TE power spectrum can break this degeneracy) . The 
projected sensitivity of 21cm emission to the IGM ther- 
mal history has been exploited to show that the 21cm 
signal can provide much more stringent constraints on a 
population of long-lived, decaying particles than is avail- 
able with CMB data 3- 

Primordial black holes (PBHs) are an excellent exotic 
physics target for 21cm observations because the mecha- 
nisms by which they heat the IGM are under tight the- 
oretical control. In this work, we compute the effects 
of a population of PBHs with masses between 10^° and 
10^5 kg on 21cm observables, including the sky-averaged 
brightness temperature and the ffuctuation power spec- 
trum. We also obtain predictions for the relic photon 
density from the population of PBHs. We take into ac- 
count the details of photon and e"^ emission by PBHs and 
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the energy-dependent absorption properties of the IGM. 
We find that future 21cm observations can provide better 
constraints on PBHs than are currently available. The 
techniques and computer code we develop can be easily 
modified to study other exotic physics scenarios in the 
Dark Ages, such as dark matter annihilation and decay, 
if the precise energy spectrum of photons and pairs 
produced by the exotic physics mechanism is known. 

The effect of PBHs on the IGM is strongly dependent 
on the PBH mass. Very massive PBHs would accrete 
matter in accretion disks and thereby emit x-rays; this 
scenario is explored in §|. We focus our present discus- 
sion on low-mass PBHs that would primarily affect the 
IGM through the Hawking radiation 0, Q of light par- 
ticles: primarily photons and pairs, and heavier par- 
ticles (and their decay products) for the smallest holes. 
The total power emission by black holes through these 
channels has been known for some time as a function of 
the black hole mass Q. Due to a coincidence between 
the Hawking temperature of PBHs that evaporate dur- 
ing the Dark Ages and a window of low optical thickness 
of the IGM to photon absorption, it is vitally important 
to know the spectrum of photons produced by the PBHs. 
The radiated photons have energies such that the proba- 
bility of absorption by the IGM, and thus their effective- 
ness in heating it, is strongly dependent on the photon 
wavelength. The emission rates as a function of photon 
energy are given by the black hole "graybody factors," 
which have been calculated. The graybody factors arise 
because black holes should not radiate with a blackbody 
spectrum, even though they act as warm bodies with a 
specific temperature. In the end, all that is needed to 
completely characterize the PBH population is the total 
density parameter ripBH and the mass of the individual 
black holes Mpbh- By comparison, while in principle 
it is possible to calculate the spectrum of photons (and 
other particles) produced in dark matter decay or anni- 
hilation, in practice this dependence is subsumed in an 
overall energy deposition rate per baryon e. There is then 
a model-dependent conversion between e and quantities 
such as the density parameter and lifetime of the relevant 
particle species. For PBHs, the energy deposition rate is 
uniquely determined once the mass of the black holes is 
specified. 

We obtain our predictions as follows. We numerically 
integrate the equations governing IGM ionization and 
temperature, including exotic sources of ener gy i njection. 
We use a modified version of the RECFAST [l3| code at 
high redshift for increased accuracy, and a simpler set of 
IGM evolution equations at lower redshifts. The energy 
injection itself is computed by using the graybody factors 
for Hawking emission and the optical depths for a vari- 
ety of IGM photon and absorption processes. Some 
of these processes redistribute photons by energy, and 
so we track the full photon population as a function of 
energy through cosmic history. We use this information 
to compute the IGM temperature and ionization history, 
and use standard techniques to compute the 21cm bright- 



ness temperature and the fluctuations. In the end, we 
obtain the complete temperature and ionization history 
of the IGM, the power spectrum and sky-averaged sig- 
nal of 21cm brightness temperature fluctuations, and a 
relic high-energy photon population. In order to obtain a 
constraint on PBHs, we compare the power spectra from 
several models to the power spectrum with no PBHs in- 
cluded and use an estimate of the measurement error for 
a realistic future experiment. 

The PBH population is currently most tightly con- 
strained for PBH masses near 

10^2 by the diffuse 
7-ray background as measured by EGRET [ll|, [l2| . The 
highest-energy EGRET photons cannot be of primordial 
origin {z > 10'^) because the IGM is not transparent at 
the corresponding energies. These photons are thought 
to originate from outside the Galaxy, but their produc- 
tion mechanism is as yet unknown. One possibility is 
that these high-energy photons arise from processes at 
redshifts lower than that of reionization, in which case it 
is unlikely that 21cm observations will be a useful probe. 
Another possibility is that these photons arise from some 
energetic process occurring during the Dark Ages. In this 
case the details of photon absorption and heating by the 
IGM will be important. If these energetic processes heat 
the IGM, then 21cm observations could provide a compli- 
mentary measurement and a consistency check on exotic 
physics mechanisms which purport to contribute to the 
high-energy diffuse 7-ray background. 

This paper is organized as follows. We review the 
physics of PBH formation and energy emission in Sec- 
tion [TTl In Section HITl we describe our computation of 
the energy injection rates from PBHs and how this in- 
fluences the mean IGM evolution, the 21cm brightness 
temperature, and the 21cm power spectrum. We give an 
overview of our results in Section IIVI We include pre- 
dictions for the ionization history, 21cm observables, and 
relic photon population for a selection of PBH models. 
In Section |V] we present the constraints on PBH pop- 
ulations that can be obtained with an ambitious 21cm 
experiment. We discuss observational prospects in Sec- 
tion |VT] and conclude in Section IVlTl 



II. PRIMORDIAL BLACK HOLE EMISSION 
A. Formation mechanisms 

Several possible formation mechanisms for PBHs have 
been discussed in the literature, and the resulting mass 
spectra are highly model-dependent. We briefly mention 
some possible PBH formation mechanisms here. For a 
more in-depth review, see, e.g., (T3| . 

As first proposed, PBHs result from fluctuations in the 
high density primordial perturbations that collapse upon 
horizon entry p^ . [l5| . This scenario is expected to pro- 
duce an extended PBH mass spectrum [lg|. As it gen- 
erally requires a spectral index > 1, this scenario is 
disfavored by the most recent observations which indi- 
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cate Us <1 (e.g., the WMAP five -year data set [l7l. Il8| ). 
Other PBH formation scenarios invoke a phase transi- 
tion or a period of a softening equation of state during 
which fluctuations would be more likely to collapse (e.g., 
[T9|). These models predict a very narrow mass spec- 
trum, as PBH formation would occur only at the time 
of the transition. Other formation scenarios invoke more 
exotic mechanisms, such as the collapse of cosmic string 
loops 0, iU, m, [2I m and of domain walls [H, HI], 
and predict a broader mass spectrum. 

In this work, we assume all of the PBHs have the 
same mass, amounting to a delta-function mass spec- 
trum. This approximation is good for some PBH for- 
mation scenarios and less so for others. It is the most 
general in the sense that its results may be easily applied 
to other models. Our results show that the constraints 
on the PBH population depend strongly on mass, but 
in a simple way. By assuming a single PBH mass, we 
can determine the mass ranges most tightly constrained 
by 21cm observations. We expect the constraints on a 
more general PBH mass spectrum to be determined, to 
a rough approximation, by the number density of PBHs 
in the mass range to which 21cm observations are most 
sensitive. Therefore assuming a single PBH mass gives a 
good indication of the constraints on more general PBH 
mass distributions. 



B. Hawking radiation 

The primary means by which black holes with masses 
< 10^'' kg would have influenced the IGM is through 
Hawking radiation. It has been known for some time that 
a Schwarzschild (uncharged, non-rotating) black hole of 
mass MpBH should radiate 0, [1| as would a warm body 
at the temperature Th, where 

^ _ hc^ _ 1.2 X 1013 K 
" " SvrGA/pBHfcs " (Mpbh/1010 kg) 
^ 1.1 GeV/fcB 

(AfpBH/lOio kg) ■ ^ ^ 

The black hole should emit all massless and nearly- 
massless particles (gravitons, photons, neutrinos), as well 
as those massive particles whose mass is substantially be- 
low ksTH- The energy emitted in neutrinos and gravi- 
tons is essentially lost, for these particles interact very 
weakly and do not affect the IGM. On the other hand, 
the IGM is affected by the energy emitted in photons and 
by the pairs emitted by smaller black holes. 

While black holes are expected to have a temperature, 
they should not emit radiation with a blackbody spec- 
trum. Instead, Hawking's calculation shows that the 
number of particles emitted with angular frequency uj 
(measured at infinity), spin s, polarization p, and angu- 
lar momentum quantum numbers Im is 

{Nsplrn{(^)) = 7T ^P'"^*- ) (2) 

exp(fitj/fcsT//) - (-1)^" 



where the Tspimi^) are graybody factors. The Tspimi^^) 
parameterize the deviation of the black hole emission 
spectrum from that of a blackbody. Physically, each 
^spirni'-^) is the probability that a particle in an infalling 
field mode described by splm is absorbed by the hole. 
Since a black hole with temperature Th must be in ther- 
mal equilibrium with a blackbody heat bath of the same 
temperature, it follows that these absorption probabili- 
ties must also determine the emission rate. In addition 
to the parameters splm, the Tspimi'-^) also depend on the 
mass of the particle species. Only particles with mass 
m <^ Th (in units where ks = c = 1) will be emitted at 
an appreciable rate, and for these particles the Tspim{^) 
are those of massless particles. In this work we only con- 
sider the emission of particles with m <C Th and so use 
the massless graybody factors. 

It is essential to include the graybody factors when 
studying the effect of black hole emission on the IGM. 
They determine the total power emission in particles of 
various spin, and therefore the evolution of the black 
hole mass. The rate of IGM photon absorption is very 
frequency-dependent, and so it is essential to know the 
photon spectrum accurately. We use the power emis- 
sion spectra computed in [9|, obtained by integrating 
the Press- Teukolsky equations for fields of different spin 
[27I [28| . A general feature of these graybody factors is 
that the emission of particles of high spin is suppressed, 
and the peak emission moves to higher energies. For ex- 
ample, the power emission per polarization is larger for 
a spin- 1/2 particle (such as a neutrino) than for a spin-1 
particle (such as a photon). For black holes with masses 
MpBH ^ 9.5 X 10^'^ kg, which are too large (and thus too 
cold) to emit ultrarelativistic pairs, most (~ 81%) of 
their energy emission is in neutrinos, which do not heat 
the IGM. Black holes with A/pbh < 9.5 x 10^^ kg, in con- 
trast, emit ~ 45% of their energy in pairs, compared 
to only ^ 9% in photons Black holes with masses 
MpBu ^ 4.5 x 10^1 kg will emit relativistic /i^ pairs, and 
black holes with progressively smaller masses will emit 
more and more massive particles. At sufficiently small 
A/pBH these will contribute to the photon, neutrino and 
electron emission by secondary decays and hadron jets. 
By a coincidence between IGM physics and black hole 
physics, the black holes that are evaporating near the 
present epoch emit their photons into the window where 
the IGM optical depth is strongly dependent on redshift. 



III. IGM PHYSICS 

The mean IGM evolution is determined by tracking the 
kinetic temperature of the gas Tr-, the ionization fraction 
Xi, and the spin temperature Ts- The evolution of Tk 
and Xi with redshift is given by differential equations, 
with initial conditions taken from standard cosmology 
just before recombination, and with additional terms to 
account for exotic energy injection mechanisms. The spin 
temperature Ts is then determined algebraically. Neutral 
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hydrogen is visible in 21cm whenever the spin tempera- 
ture Ts differs from the CMB temperature, Tcmb, and 
due to the contrast, the signal will show up in either emis- 
sion or absorption. The sky-averaged brightness temper- 
ature Tf, therefore carries information about the mean gas 
temperature and ionization state of the universe. How- 
ever, because of the many bright foregrounds expected for 
21cm experiments (discussed in more detail in § VI A|) . an 
absolute all-sky signal will be difficult to detect. An eas- 
ier target is a statistical detection of the power spectrum 
of the 21cm brightness temperature perturbations. By 
comparing the power spectrum of the 21cm signal with 
and without a contribution from PBHs, we can retain 
the ability to discriminate among models within the lim- 
itations of realistic experiments. We discuss the power 
spectrum calculation in ijlllBI 



A. Mean IGM evolution 

The homogeneous IGM is described by the mean ki- 
netic temperature Tk, the spin temperature Ts, and the 
ionization fraction Xi. The mean kinetic temperature of 
the gas evolves according to [1, [1^ 



dT, 



K 



dt 



= -2H{z)Tk 



(Tcmb — Tk) 



(3) 



where H{z) is the Hubble parameter at redshift z, rji — 
1 + fHe + Xi with .fee the helium fraction (defined by 
riHelViH + nHe]) [HI, and 



^ 3meC 
^ ^cttUcmb 



(4) 



with rrie the electron mass, ctt the Thomson cross section, 
and C/cMB the CMB energy density at redshift z. The 
first term on the right-hand side of Equation ^ accounts 
for the redshifting of kinetic energy with the expansion of 
the universe. The second term includes the effect of IGM 
heating by scattering of CMB photons from hydrogen 
ions. The final term in ([3]) takes account of heating by 
exotic energy injection into the IGM. The parameter e 
is an energy injection rate per baryon, and Xh is the 
fraction of the energy that goes into heating the IGM. 
We describe in more detail below how e is determined for 
the PBH population. 

The ionized fraction Xi obeys [2^ 



dXj 

"dT 



-a[TK)x1nH 



(5) 



where the first term on the right-hand side includes 
hydrogen recombination through an effective coefficient 
a{TK), which depends on the kinetic temperature of the 
gas. Following we use the Case B coefficient from Ta- 
ble I of [13] ■ The second term in Equation ([5]) includes 
the ionizations produced by an exotic energy injection 
mechanism. As in Equation ([3]), e is the energy injection 



rate per baryon, Xi gives the fraction of energy going into 
ionizations, and Eion = 13.6 eV is the ionization thresh- 
old of hydrogen. 

The spin temperature of the neutral hydrogen gas de- 
pends on competition among coUisional excitation [3l| . 
heating by CM B p hotons, and interactions with Ly-a 
photons [H,!!!,!!^. The temperature itself is defined by 
the ratio between the occupation numbers ni and tlq of 
the singlet and triplet hyperfine spin states, through 



"1 Q ( 

— = 3exp - — 
no V 



(6) 



where = 0.068 K is the energy splitting between hy- 
perfine states in temperature units. In equilibrium, the 
spin temperature Ts is given by 



1 



1 



Ts 



J CMB J^K J c 



(7) 



where the parameters Xc and Xa are the collisional and 
Wouthuysen-Field [H, [111 coupling parameters, respec- 
tively, and Tc is the effective color temperature of the 
radiation field [s^, [1^ • We omit the Wouthuysen-Field 
coupling in this work, since we only expect it to be impor- 
tant once the first luminous sources turn on at compar- 
atively low redshift, at which point their heating effects 
will easily swamp the contribution from PBHs. The col- 
lisional coupling coefficient Xc is 



4Ki_o(rg)ngr^ 
SAiqTcmb 



(8) 



where Aiq = 2.87 x 10""'^^ s~^ is the spontaneous decay 
rate of the hyperfine transition. The parameter ki^q 
describes spin excitation through atomic collisions; its 
value has been tabulated as a function of gas temperature 
HE [13. For the range 1 K < IV < 300 K, we use the 
values in column 4 of Table II in [13 . For Tk > 300 K, 
we use the fitting formulae suggested in [13 for the data 
presented in (36j . 

Using Equations (O and Q to evolve Tk, Xi and 
Ts, we calculate the differential 21cm brightness relative 
to the background [J, [11] : 



28 mK 



-1/2 



CMB 



(9) 



(10) 



0.033/ V0.27 

/ Ts-T, 
10 7 ^ Ts 

where we assume r <C 1 in the second line. 



B. 21cm power spectrum 

Measurements of the power spectrum of fluctuations in 
the 21cm signal will make it possible to extract statistical 
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information about the structure in neutral hydrogen even 
if detailed imaging or tomography remain beyond the 
reach of experiments. 

As evident in Equations (O and ^ above, the bright- 
ness of the 21cm signal depends on the density, temper- 
ature and ionization state of the hydrogen gas. To trace 
the shape of the power spectrum as a function of redshift, 
we must track not only the density perturbations in the 
baryonic matter, but also the fractional perturbations in 
the ionization state and the temperature of the hydrogen 
gas and the factors that influence the couplings among 
them. 

The fractional perturbation to the brightness temper- 
ature is defined as 



(11) 



The Fourier transform (52i(k) can be written as a sum 
of each contribution weighted by expansion coefficients 
related to the various couplings [1, 0] : 

52i(k) = ip + + finSH + Mo. + Mt, (12) 

where the 5i are fluctuations in overdensity ((5), neutral 
fraction {5h), Lyman-a coupling strength {5a) and tem- 
perature {5t)- The expansion coefRcients are 



/3 

Pa 



Xtot{l + Xtot) 



„HH 



X 



eH 



Xtot {I + Xtot) 



Xtot(l + Xtot) 
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1 



eH 



X4ot(l + Xtot) 



dluKf^ „„dlnK^," 



(13) 

(14) 
(15) 

(16) 



where the Xc are collisional coupling coefRcients [39j . 
for which Xc = x'^^^ -\- x^^ and Xtot = Xc + Xa- 
The parameters k^q and are rate coefficients for 

spin de-excitation from electron/hydrogen and hydro- 
gen/hydrogen collisions, respectively 01 ■ 

The power spectrum of 21cm fluctuations contains con- 
tributions proportional to the mean brightness temper- 
ature of the 21cm line over the whole sky as well as 
to the fluctuations in the spin temperature j40| . Tak- 
ing a simplified estimate in which we consider only the 
growing-mode perturbation on large scales, we can write 
all perturbations as pure time-delay perturbations, which 
means each separate perturbation, and hence the 21cm 
power spectrum, is proportional to the matter power 
spectrum Pss- As we discuss in more detail below, this 
approximation leads to a small error when the brightness 
temperature is near zero, but provides sufficient accuracy 
for our purposes. Then, the 21cm power spectrum can 
be written 



P2i{k,ti) = 5T^{(3' + lx')Pss{k), 



where /i = cos(0) accounts for the angle between the 
wavevector k and the line of sight, and we have 



13' = (3 + /?T5T 



(1 - Xi) 



(18) 



Here, gi{z) = Si/S and gxiz) = 6t/S are defined for 
convenience to compare the fiuctuations in the ioniza- 
tion and temperature to the matter overdensity, and 9u 
encodes the uniformity of the energy deposition. It is de- 
fined such that the energy deposition rate is proportional 
to {l + 9uS) so that 9u = for uniform deposition. Recall 
that we are neglecting the contribution of Wouthuysen- 
Field coupling in this work, so we are effectively imposing 

/?a =0. 

In practice, the power spectrum is calculated by evolv- 
ing gi and gx with redshift along with the mean global 
properties of the IGM. The evolution equations are Q 



dffT _ .gT-2/3 Xj grTj - gi{Tj - Tk) 
dz ^ 1 + z T]it-y Tk{1 + z)H{z) 

2 T]2e {l-0u)+gT 
-Xh- 



dgi 
dz 



iTliksTK {l + z)H{z) 
gi aXiTiHil + gi + a'^r)) 



l + z {l + z)H{z) 

e (1 - Qu) + gi 



Eion Xi (1 -I- z)H{z) 



(19) 



(20) 



(17) 



where we have included bars over quantities that are 
global averages, for clarity. The quantity a is the hy- 
drogen recombination coefficient and a' = dlna/dlnT^ 
i. 



C. Energy injection and the photon spectrum 

The PBHs influence the IGM through the direct emis- 
sion of photons and pairs. Given the black hole emis- 
sion power in these channels we must take into account 
the details of photon and absorption by the IGM in 
order to give the power absorption per baryon e which 
appears in the evolution equations for Tk and Xi. 

There are a number of processes that allow photons to 
deposit energy in the IGM, all of which have frequency 
dependencies, which themselves change with redshift and 
with the ionized fraction. We follow the tabulation of op- 
tical depths given in [4l[ for the following processes: pho- 
toionization, Compton scattering, pair production from 
atoms and ions, photon-photon scattering, and single and 
double pair production from CMB photons. Individual 
rates for each process are shown at a fixed redshift in 
Figure m and the total photon-IGM interaction rate as a 
function of redshift is shown in Figure^ At all redshifts, 
below roughly 100 eV, the optical depth is dominated by 
photoionization, while at high energies it is dominated by 
pair production from CMB photons. Between the two, 
the absorption of photons depends strongly on energy. 
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Photon absorption at z=300 Photon transparency window 




proper photon energy (eV) 



FIG. 1: Optical depths per time for various photon-IGM pro- 
cesses, in units of the Hubble time, at z — 300, assuming a 
neutral IGM. These include processes which deposit energy 
directly into the IGM (pair production and photoionization) , 
processes which redistribute photons (27 27) and ones that 
do both (Compton). At very low energies, photoionization is 
the dominant process; at very high energies, pair produc- 
tion dominates. 



At lower redshifts, a "transparency window" for photon- 
IGM interactions emerges [a,|H,|42l. The "bump" visible 
in the plot is due to Compton scattering, and photons 
that fall into the trough at slightly higher energies are 
largely unabsorbed by the IGM. 

To enable us to fully account for the IGM/photon 
physics, we follow the photon population in 10^ co- 
moving energy bins between 1 and 10^^ eV. (We bin by 
comoving energy so that redshifting alone does not move 
photons from bin to bin.) Careful accounting is nec- 
essary because of a numerical coincidence between the 
transparency window (the trough between 10^ — 10^ eV 
at z = 300) and the Hawking temperature of long-lived 
primordial black holes. The black hole mass that is most 
interesting from the perspective of IGM physics is that 
corresponding to black holes that are just evaporating 
away at the present epoch, for these black holes deposit 
almost all of their initial mass-energy in the IGM. A black 
hole that is just evaporating today had an initial mass of 
5 X 10^^ kg 0, corresponding by Equation ([1]) to a tem- 
perature of 20 MeV, which is nearly constant until the 
very end of the hole's life. Referring to Figure [5] reveals 
that these black holes are emitting a substantial fraction 
of their photons near the transparency window. There- 
fore, most photons emitted by these black holes will not 
be absorbed for many Hubble times, if at all. In this case. 



FIG. 2: Total optical depth for photon-IGM processes at var- 
ious redshifts. For photon energies where (dr/di)/// < 1, the 
IGM is essentially transparent. The Hawking temperature 
Th for various PBH masses is also shown. Because of the 
graybody factors, the peak in energy emission is at ener- 
gies slightly above ksTn, and the peak in the photon energy 
emission higher still. 

treating black hole emission as absorbed "on the spot" 
would be a poor approximation. To properly treat black 
hole emission near this mass range, then, we follow the 
full photon population over time. 

Different photon/IGM processes deposit energy in the 
IGM in different ways. The pair production processes 
(from atoms or ions, and single and double pair pro- 
duction from CMB photons) all deposit their energy by 
converting it into pairs. We count all of the energy 
going into the pair (including the rest masses of the 
themselves) as injection into the IGM via the e pa- 
rameter above. When deposited in the IGM, the energy 
is partitioned between heating, ionization, and HI exci- 
tation. The fractions of the initial electron energy that 
go into heating and ionization are denoted by Xh and 
Xi, which were introduced in ^ and (O. These frac- 
tions have been computed in detail [isj ; for our work we 
employ the approximations [1, 112] 



We treat photoionization within the same framework. 
The initial photoionization event is different than heating 
or ionization caused by the energetic electrons produced 
by pair production. Nonetheless, except for photons near 
the photoionization threshold of 13.6 eV, the majority of 
the initial photon energy is given to the secondary pho- 
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toelectron, which then heats and ionizes subsequent hy- 
drogen atoms and ions in the same way as the electrons 
produced by pair production. Tracking the full photon 
population allows us to properly treat photons whose first 
interaction occurs long after they are created, but for all 
of the pair production processes, we assume that the en- 
tire photon energy is lost on this first interaction. 

In contrast to the pair production processes, other 
processes largely redistribute photons within the photon 
population. These photons do not deposit their energy 
in the IGM directly, though they may do so after down- 
shifting to lower energies. The main example of such 
a process is two-photon scattering 77cmb 77 of ener- 
getic photons from the CMB. We assume that each inter- 
action evenly divides the initial energy between the two 
photons, and so this produces a cascade of photons at 
successively lower energies. These may eventually inter- 
act, or fall into the transparency window and be lost. To 
a certain extent, Compton scattering also redistributes 
photons among the energy bins. Unlike the other pro- 
cesses we consider, each Compton scatter may only re- 
duce the energy of the photon by a small amount, and so 
we cannot make the approximation that all of the initial 
photon energy is deposited in the IGM. We deposit an 
energy-dependent fraction (given in [4l|) of the photon's 
energy in the form of hot electrons, and decrement the 
energy of the photon itself. Thus Compton scattering 
both heats and ionizes the IGM and modifies the photon 
population. 



IV. NUMERICAL RESULTS 

We have developed a modular code that simultane- 
ously: (a) evolves the PBH population; (b) calculates 
the energy output in photons and other particles; (c) ac- 
counts for partial absorption and free streaming due to 
the transparency window (evolved with redshift); (d) par- 
titions the energy injection into heating and ionization; 
(e) evolves the sky-averaged quantities; and (f) calculates 
the perturbations in heating and ionization for the 21cm 
power spectrum. 

For initial conditions, we use a modified version of the 
RECFAST tlfl] code starting at redshift z = 10'', and 
then use that code to continue to follow the evolution of 
the mean temperature and ionization and their pertur- 
bations up to a redshift z — 300. After that point, we use 
a simplified IGM evolution code to implement Equations 
® , ®, HH), and (HOI)- Both RECFAST and the simpler 
IGM evolution code use an additional code module which 
tracks the photon population, and which computes the 
energy injection rate from photons and pairs which 
feeds back into the IGM evolution equations. 

For the purpose of this calculation, we begin with a 
population of PBHs at a single mass and a given num- 
ber density. The possible spectra of masses of PBHs is 
discussed in i jll Al 

We follow the evolution of the population of PBHs as 



they lose mass to Hawking radiation, injecting their en- 
ergy into the IGM. Based on the PBHs' effect on the 
21cm power spectrum in a range of redshifts and the ex- 
pected observational limitations of future experiments, 
we are able to obtain a limit on the density of PBHs at 
each initial mass in our range of interest. 

For most of the mass range of interest for 21cm ex- 
periments, the PBHs have evaporated completely before 
the present day. As can be seen in Equation ([1]), the 
Hawking temperature diverges as the PBH mass goes to 
zero. There has been much discussion in the literature 
as to the phenomenology of the final stages of Hawking 
evaporation, with some suggesting that the PBHs may 
leave a Planck- mass relic [4lli45l. |46|. and others propos- 
ing a connection between the late stages of PBH evapo- 
ration and signatures of extra dimensions For our 
purposes, the details of the final moments are not impor- 
tant, since only a very small fraction of the PBH's total 
energy is emitted during that time. We "turn off" the 
Hawking evaporation when the PBHs are within one year 
of complete disappearance. 

We find there are several regimes of PBH mass ranges 
that have distinct effects on the evolution of the IGM: 

• -MpBH ^ 5 X 10^° kg: In this regime, the evapo- 
ration occurs early enough to alter the late stages of 
recombination, but these PBHs are not well suited 
to constraints with Dark Age 21cm observations. 

• 5 X 10^° kg < MpBH < 10^^ kg: PBHs in 
this range of masses evaporate in a range of red- 
shifts (30 S ^ S 90) during which the fiducial 
sky-averaged 21cm brightness temperature signal 
is visible in absorption against the CMB. At suf- 
ficient number densities, these PBHs would raise 
the brightness temperature in 21cm, effectively de- 
creasing the strength of the 21cm signal during this 
time. An anomalously small (in absolute value) 
21cm brightness temperature on the sky would in- 
dicate extra energy injection in this regime. 

• AipBH ^ lO-'--'- kg: The evaporation of PBHs with 
initial masses near 10"'^^ kg have the most dramatic 
effect on the expected 21cm signal. The peak in 
their energy injection would occur at sufficiently 
low redshifts that, at high densities, they could 
raise the spin temperature above the temperature 
of the CMB, making the 21cm signal appear in 
emission rather than in absorption. For the small- 
est masses, at number densities not yet ruled out 
by current data, the emission signal can be several 
mK. If the all-sky signal were to switch to emission 
before the effects of star formation were expected 
to significantly alter the signal (stars can bring the 
21cm signal into emission for redshifts less than 
z ~ 15 [48|), this would be a strong indication of 
exotic physics in the Dark Ages Q . 

• 10" kg < MpBH < 10" kg: PBHs created 
in this mass range would still exist today, as their 
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evaporation timescale is longer than the current age 
of the universe. For these PBHs, the Hawking radi- 
ation can be well approximated by a constant rate 
of energy injection; their effect on the IGM would 
be very similar to that of decaying dark matter (see 
for a discussion of the effects of such a scenario) . 
Since the power and Hawking temperature are both 
small for such massive PBHs, their effect on the 
IGM would be less pronounced than that of their 
less massive counterparts. They would also lack 
the sudden, short-lived increase in signal that can 
be seen for PBHs that complete their evaporation 
in the Dark Ages. 

• MpBH > 10^^ kg: At very high masses, the 
Hawking temperature is too low to allow for the 
production of pairs. These black holes emit 
cooler photons which do not efficiently inject energy 
into the IGM. Therefore a greater number density 
of these black holes can be tolerated without chang- 
ing the mean IGM evolution significantly. 

For black holes with masses above 10^^ kg, which do 
not emit pairs, the energy injection rate is roughly 
constant. This allows for a direct comparison with de- 
caying dark matter models in which the dark matter life- 
time is much longer than the age of the universe, which 
give a roughly constant energy injection rate per baryon. 
Therefore we can find a dark matter decay model corre- 
sponding to a given large-mass PBH evaporation model. 
For cases in which all the dark matter is in the form of 
decaying dark matter particles, the correspondence is: 

/ MpBH \ / npBH \ 

V 1015 kg) \ 1.2 X 10-42 m-3 J 

= (3.18x10-30 s-i) ^^^^ 

for MpBH > lO^"' kg, where Fx is the decay rate of 
the dark matter particle. As a corollary, observational 
consequences in 21cm of PBH populations are a function 
of Mpi^QpBH for MpBH > 10^"* kg. For smaller PBHs, 
the emission of pairs dominates the energy injection, 
and this scaling no longer applies. 

To determine the detectability of PBHs evaporating in 
the Dark Ages, we simulate their effects in three poten- 
tially observable quantities. The first is the overall ion- 
ization state of the IGM, the second is the high-redshift 
21cm all-sky brightness temperature, and the third (and 
most relevant for near-term experiments) is the 21cm 
power spectrum. We discuss each in turn in the following 
subsections. 

To illustrate the effects of PBHs in different mass 
regimes, we will use three example PBH models, whose 
effects on the IGM are plotted in Figures El HI and [S] The 
parameters of the three cases are described in Table HI 



A. Mean ionized fraction 

As the PBHs evaporate, the energy they inject into the 
IGM will increase the hydrogen ionized fraction. To ex- 
amine the importance of this effect, we simulated the ion- 
ization history of the IGM, following Equation ^ with 
the energy deposition for ionization determined by Equa- 
tion (PT|) . For the models we considered, which would af- 
fect the 21cm power spectrum at potentially detectable 
levels (see HIV C[) . the ionization fraction did not increase 
to a level sufficient to reionize the universe; generally, 
ionization fractions of the order of 10"^ are only reached 
at redshifts low enough that the formation of stars and 
galaxies would likely start to become important. We plot 
some example ionization histories in Figure [3l 

The models we plot here are chosen to reflect the range 
of phenomenology that can be expected from PBH evap- 
oration in the regimes in which it would strongly affect 
the IGM evolution. In addition to the no-PBH case, we 
plot two cases in which the PBH evaporation completes 
during the Dark Ages, giving the strongest signal as the 
temperature and power of the Hawking radiation increase 
greatly in the final stages of evaporation, and one case in 
which the evaporation is not complete by the present day. 
The parameters of these models are summarized in Table 
|I1 In Figure [31 the peaks in each of the two lowest-mass 
models occur at the redshift at which the evaporation 
finishes in each model. 

One may ask if the alteration in the ionization his- 
tory would be detectable as a change in the integrated 
Thomson scattering optical depth measurement (tlss) in 
GMB experiments. For the models we focused on, with 
potentially detectable differences in the 21cm power spec- 
trum (as determined by equation [23]), we find that the 
change in tlss is at best on the order of 10"^^ whereas 
the contribution to tlss from the reionization of the uni- 
verse is two orders of magnitude higher. Current exper- 
iments cannot reach this precision, and uncertainties in 
the process of cosmological reionization would likely pre- 
vent exotic physics from being constrained in this way 
in the future: for example, a change in the reionization 
redshift from 6.1 to 6.0 would result in a change in tlss 
of order 10"^^ and non-instantaneous reionization sce- 
narios complicate the picture further. The fact that the 
optical depth measurement is an integrated signal over 
the history of the universe since recombination gives it a 
disadvantage over future 21cm experiments, where mea- 
surements will be possible at a large number of redshifts. 



TABLE I: Example primordial black hole models for Figures 
HHandEl 



Model 


MpBH (kg) 




Plot key 


A 




1.3 X 10"^ 


green, dot-dashed 


B 


10" 


8.2 X 10"!^ 


blue, dashed 


C 


5 X 10^° 


2.7 X 10^1^ 


magenta, dotted 
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FIG. 3: The mean ionized fraction vs. redshift for the fiducial 
case of no PBH energy injection (thick sohd line) and three 
PBH cases, described in Table H] A/pbh = 10^^ kg at Upbh = 
1.3 X 10"'^ (Model A; green, dot-dashed), Mpbh = 10 kg 
at fipBH = 8.2 X 10"^^ (Model B; blue, dotted) and Mpbh = 
5x 10^° kg at fipBH = 2.7x 10"^^ (Model C; magenta, dotted). 
These models correspond to marginally detectable deviations 
in the power spectrum, as discussed in Equation (I23|l in i]IV CI 



With CMB polarization measurements from the Planck 
satellite, we expect to be able to constrain energy in- 
jection rates of ~ IQ-^^ eV/s/baryon at z = 1000 
whereas with future 21cm power spectrum measurements 
we can expect to do around three orders of magnitude 
better 0], and to measure at several redshifts. 

Since the 21cm signal depends on the ionization state 
and the temperature of the gas, the features in the ion- 
ization history are echoed in the 21cm brightness tem- 
perature signal, which we discuss next. 



B. Sky-averaged 21cm brightness temperature 

The ultimate realization of the power of 21cm obser- 
vations would be in the form of 21cm tomography of the 
high-redshift intergalactic medium. A radio telescope ob- 
serving the full sky at 21cm(l-|-z) could, in principle, have 
a complete map of the neutral hydrogen at that redshift, 
and could obtain a 3D image of structure formation by 
combining slices at each z. While there are significant ob- 
servational challenges involved in true tomography (some 
of which we discuss in WII) . even a sky-averaged mean 
signal of the 21cm brightness temperature, which would 
contain information about the overall state of the neutral 
hydrogen at a range of redshifts, would significantly im- 
prove upon current probes of the IGM at high redshifts. 

Following Equation we track the mean 21cm 

brightness temperature as a function of redshift in our 
simulations. Although there are still many observational 



FIG. 4: The sky-averaged brightness temperature vs. redshift 
for the fiducial case of no PBH energy injection (thick solid 
line) and three PBH cases from Table IT] as in Figure [3] above. 
The thin horizontal line is at a 21cm brightness temperature 
of zero (no signal) . 



difficulties for this measurement (including major fore- 
grounds and challenges to instrument design, discussed 
in more detail in WI|) . this result gives an impression of 
the effects of PBH energy injection on the 21cm observ- 
ables, and experiments currently under way may be sen- 
sitive to strong gradients in the sky-averaged signal. In 
Figure m we plot the brightness temperature induced by 
the no-PBH case as well as by the three example models 
plotted previously in Figure [31 

The two low-mass cases, in which the final evapora- 
tion occurs in the Dark Ages, show conspicuous spikes 
in the brightness temperature at the redshifts of evapo- 
ration. For the lowest-mass case, Mpbh = 5 x lO^*' kg 
(Model A), the effect is to decrease the 21cm signal in 
absolute value by tens of mK; this discrepancy should be 
noticeable if future experiments can measure the mean 
signal to sufficient accuracy. The signal from the inter- 
mediate case (Model B) would potentially be much more 
striking, as the signal would switch from absorption to 
strong emission at a redshift above that at which ordi- 
nary star-formation is expected to bring the spin tem- 
perature above the CMB temperature. The highest-mass 
case (Model C), in which evaporation has not occurred by 
the present day, would also significantly alter the mean 
signal over a range of redshifts, though there would not 
be as dramatic a discrepancy. 

It must be noted that the "spikes" in the ionization 
and brightness temperature signals depend upon the si- 
multaneous evaporation of all PBHs — this kind of signal 
would only occur if the PBHs were formed with a very 
narrow mass spectrum. For a broader range of masses, 
the peaks would be smoothed out and the signals less 
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FIG. 5: The 21cm power spectrum for the fiducial case of no 
PBH energy injection (thick solid line) and the three PBH 
cases from Table U] as in Figures [3] and |4] above. 



striking. Therefore the absence of a spike in the sky tem- 
perature would not necessarily disfavor the existence of 
PBHs with a broad mass distribution. 



C. 21cm power spectrum 

Given the observational difficulties inherent in 21cm 
tomographic imaging of early structures and in all-sky 
brightness temperature measurements, some of the near- 
term experiments are instead aimed at detecting the 
21cm fluctuations statistically by measuring the 21cm 
power spectrum. Realistic near-term experiments should 
have the capability to use the 21cm power spectrum to 
study the process of reionization [i^ and estimate cos- 
mological parameters [s^l, and we find that PBH evap- 
oration can be strongly constrained by more ambitious 
future radio interferometers. 

The ability to distinguish models including exotic 
physics from the standard IGM evolution depends upon 
how precisely the power spectrum can be determined. 
We use an estimate of the error in the power spectrum 
measurement for a hypothetical future experiment to test 
the observability of the alteration in the IGM evolution 
due to our PBH models. The error depends on several 
different parameters of the experiment in question. An 
ideal radio interferometer would have a long maximum 
baseline (i?max), a high effective area {A^ = N x A, 
where N is the number of receivers and A is the area of 
each), and a high covering factor {fcov = ^e/(7ri?niax))- 
A large bandwidth (_B) and a long integration time (tint) 
also decrease the error. The frequency of the observa- 
tion comes into the equation in two ways, and in both 
cases, higher redshifts give higher errors. The error is 
proportional to (1 -I- z), but also to the sky temperature 



Tsky, which is a function of frequency. The foreground is 
dominated by Galactic synchrotron, which roughly fol- 
lows Tsynch ^ (i^/200 MHz)"2-^ K for the frequencies of 
interest. 

For the purpose of placing a constraint, we consider an 
experiment with parameters similar to estimates for the 
proposed Square Kilometer Array (SKA), but operating 
at a much lower frequency in order to probe the higher- 
redshift universe 5l| . We compare the parameters of the 
considered low-frequency experiment with those of SKA 
in Table im We assume an integration time of 1000 hr for 
each, and observations at wavenumber k = 0.04 Mpc^^. 



TABLE II: Experimental parameters for proposed future low- 
frequency observation compared with expected SKA param- 
eters. 





SKA 


Future experiment 


-^max 


5 km 


5 km 


fcov 


0.01 


0.25 


bandwidth 


2 GHz" 


50 MHz 


min. frequency 


100 MHz 


29 MHz 


max. 1 + z 


14 


50 


Tsky 


660 K 


10* K 



"The bandwidth used for an individual observation is frequency- 
dependent; see |5ll for details. 

From the error at a wavenumber k is given by: 



fc3(5P2i 0.1 mK /0.25\ /5 km 



27r2 



k 

0.04 Mpc 
/lOOO hr^ 



GOV / \ ^max 

3/4 



( ^sky 



VIO^ K 

50 



50 MHz\ 



1/4 



(23) 



The data are assumed to be binned in segments of loga- 
rithmic length ek. 

In Figure [51 we plot example power spectra from our 
PBH models alongside the fiducial no-PBH model to il- 
lustrate how the power spectra change with PBH energy 
injection. The zeros in the power spectra plots corre- 
spond to the 21cm signal changing from absorption to 
emission, since, from Equation (jl7p . for our approxima- 
tion, the power spectrum is proportional to the square of 
the 21cm brightness temperature. 

We have made some approximations in order to com- 
pute these power spectra, and here we can comment on 
their accuracy. Essentially we have taken the 21cm power 
spectrum to be proportional to the mean 21cm brightness 
temperature Tb- However there is also a contribution 
from fluctuations in the spin temperature Tg. Normally 
these are small but can become important when the 21cm 
mean signal switches from absorption to emission and 
passes through zero. We can use the more accurate calcu- 
lation of 1401 to estimate our error. The fractional error, 
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comoving photon energy (eV) 

FIG. 6: Evolution of the photon population for Model A of 
Table U with A/pbh = 10^^ kg. The solid curves, from left 
to right, are z = 4000, 1000, 300, 100 and 30. The dashed 
curve is the photon population at z = 30 with all IGM-photon 
interactions turned off, given for comparison. 
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FIG. 7: Evolution of the photon population for a model with 
MpBH = 10" kg and J7pbh = 1.3 x 10"^ The solid curves, 
from left to right, are z = 4000, 1000, 300, 100 and 30. The 
dashed curve is the photon population with all IGM-photon 
interactions turned off, given for comparison. 



relative to the no-PBH 21cm power spectrum, is approx- 
imately 



Ts At, 
Ts - Ahi 



(24) 



using the notation of [40]. This is the ratio between 
the fluctuations in the spin temperature (which we have 
neglected) to those proportional to the mean brightness 
temperature (which we have included) . Using the results 
of [40] for the standard cosmology (no PBHs) gives a 
fractional error of ~ 1/7. While not suitable for percent- 
level precision observations, our approximations can be 
trusted to indicate when PBHs lead to a substantial de- 
viation from the expected 21cm signal, and give reliable 
results at the ~ 10% level. 

Applying the criterion of Equation (^5]) to these power 
spectra, we can estimate the ability of future 21cm ob- 
servations to place limits on populations of PBHs. We 
discuss our results in the Section Ivl 



two models follow a similar evolution pattern. At very 
high rcdshift (z > 1000) the black holes emit photons 
near the Compton peak, which are subsequently down- 
scattered to lower energies. After recombination, the low- 
energy photons rapidly deposit all their energy in the 
IGM through photoionization of neutral hydrogen and 
disappear from the spectrum. As the redshift decreases 
further, most photons are emitted into the transparency 
window. The late-time photon population therefore re- 
laxes to the unabsorbed Hawking spectrum, with a tail 
at lower energies from Compton downscattering. 

The models shown in the Figures are chosen to be on 
the constraint surface for 21cm observations. For these 
models, the relic photon population would be swamped 
by conventional astrophysical processes in the relevant 
energy range. Hence 21cm observations could detect the 
presence of these PBH populations that would be unob- 
servable from diffuse photon background measurements 
alone. 



D. Relic photon population 

As a by-product of the IGM computations, we ob- 
tain predictions for the relic photon population produced 
by the PBHs. These photons have an energy spectrum 
which is the Hawking spectrum, modified by the gray- 
body factors, convolved over redshift, and then processed 
by IGM-photon interactions. 

Plots of the photon population as a function of rcdshift 
are given in Figures [H] and [T] The populations in these 



V. CONSTRAINTS ON PBH POPULATION 

The power spectrum of 21cm brightness temperature 
fluctuations will be a powerful probe of exotic physics 
in the Dark Ages, before the IGM is signiflcantly heated 
and ionized by the first stars and galaxies. We apply an 
estimate of the error in the power spectrum, from Equa- 
tion (j23|) . to our simulated power spectra for models of 
PBHs injecting energy into the IGM through Hawking ra- 
diation. We find that future experiments will be able to 
place very tight constraints on PBHs that evaporate dur- 
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MpBH (kg) 

FIG. 8: The solid line shows our expected constraints on 
PBH mass and density parameter from a 21cm observational 
program with parameters described in Equation ((23}. The 
blue dotted line shows the diffuse 7-ray background constraint 
from [53 | for comparison. PBH models with masses less than 
~ 5 X 10^° kg are difflcult to constrain with the 21cm obser- 
vations we are proposing because they evaporate before the 
Dark Ages, and may interfere with the recombination process. 



ing the Dark Ages, improving on current limits by several 
orders of magnitudes in certain regimes, and demonstrat- 
ing the feasibility of using 21cm observations to constrain 
exotic physics. 

In Figure m we plot the constraint based on the simu- 
lated power spectra, using Equation ([^5]) as an estimate 
of distinguishability, along with the PBH constraint from 
the diffuse gamma ray background [52]. 

For the mass range in which the evaporation is not 
complete before the present day, the energy injection 
from evaporation remains roughly constant during the 
Dark Ages, and the constraint on the density parame- 
ter (ripBH) hi PBHs is proportional to A^pbh- For the 
case of a constant energy injection rate, the fact that the 
Hawking radiation power is proportional to Mpjjj im- 
plies that the combination Mpg^OpBH determines the 
energy injection rate, so the scaling in that regime is as 
expected. 

For intermediate masses, where the evaporation fin- 
ishes during the Dark Ages, somewhat tighter constraints 
are possible due to the sharp increase in the power and 
Hawking temperature during the final stages of evapora- 
tion. 

At the smallest masses, the evaporation occurs be- 
fore recombination is complete; this case is not well con- 
strained with Dark Age probes. 



VI. OBSERVATIONAL PROSPECTS 

The prospect of using 21cm observations to detect the 
transition from a neutral to ionized IGM has motivated 
a surge of activity in the observational sector, resulting 
in several ambitious low-frequency radio interferometry 
projects, some of which are already taking data. 

A review of the progress of current and planned obser- 
vational projects can be found in [49[. For our purpose in 
discussing a possible future constraint on exotic physics 
in the Dark Ages, we imagine an ambitious future radio 
interferometry experiment that will measure the 21cm 
power spectrum at 0.1 mK precision at a redshift of 50. 
The parameters of such an experiment (see Table [H]) are 
chosen to demonstrate the potential of 21cm observations 
rather than to point to a specific proposed experiment; 
our imagined experiment would be somewhat more ad- 
vanced than the proposed SKA. 

Sensitive study of neutral hydrogen in the Dark Ages 
with 21cm observations will require overcoming chal- 
lenges both in instrumentation and in astrophysical fore- 
ground removal. In terms of the instrumentation, it will 
be important to have tight control over the polarization 
response of the instrument and to carefully deal with 
frequency-dependent side-lobes [49| . 

An even bigger challenge to this observation will be 
the strong astrophysical and terrestrial foregrounds. We 
discuss these briefly in the next section; for more details, 
see, e.g., 0, IH, HI Hi US HI ■ 

A. Foregrounds 

The radio sky at the frequencies of high-redshift 21cm 
observations is very bright; at around 200 MHz (corre- 
sponding to 21cm observations at z ^ 6), the temper- 
ature of the quietest parts of the sky can be hundreds 
of Kelvins, compared to an expected 21cm fluctuation 
signal on the order of tens of mK . Radio synchrotron 
emission from our Galaxy is the dominant contribution to 
this foreground; other contributions are free-free emission 
from our Galaxy and extragalactic contributions. The 
extragalactic component could be particularly difficult 
to correct for in spatial measurements because the angu- 
lar fluctuations are unknown and will swamp the 21cm 
signal [53]. 

Frequency information obtained along individual lines 
of sight may alleviate the foreground problem to some 
extent. Synchrotron and free- free emission should have 
smooth powcrlaw spectra in frequency space, whereas the 
high-redshift 21cm signal from neutral hydrogen struc- 
tures is expected to have more small-scale structure (see, 
e.g., [5^). After point sources are excised, a smooth spec- 
trum can be flt and subtracted, leaving only the 21cm 
fluctuations. This would generally be done on a pixel- 
to-pixel basis, using a fitting function for the smooth 
foreground spectrum. It would be especially effective in 
cases, such as with the SKA, where individual HII regions 
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may be imaged. These fully-ionized regions will have no 
signal in 21cm, and so will allow for a direct probe of the 
foregrounds 

Fluctuations in the ionosphere can also be problem- 
atic for low-frequency measurements, distorting images 
through refraction. These distortions will likely be cor- 
rectable except during periods of unusually high iono- 
spheric activity [i^. Below frequencies of a few MHz, 
one approaches the plasma frequency of the ionosphere, 
and it becomes opaque. This effectively prevents ground- 
based 21cm observations at redshifts higher than z ~ 70. 

One of the most difficult foregrounds to correct for will 
be terrestrial radio interference, though there is much on- 
going computational work being devoted to this problem. 
Radio, television, and mobile phones all operate in the 
desired frequency range, as it is not a protected region of 
the spectrum. For this reason, many of the current and 
upcoming observational efforts are occurring in remote 
locations with low population density. To completely de- 
feat the foregrounds from both the ionosphere and radio 
interference, some have suggested constructing observa- 
tories on the far side of the moon, with the lunar mass 
acting as a radio shield [13, 113 ■ 

Tight instrumental control and careful foreground re- 
moval will be the biggest barriers to the use of 21cm ob- 
servations for cosmology or early universe physics, and 
the ultimate goal of detailed imaging of early structure 
formation in neutral hydrogen may be thwarted by these 
challenges for the foreseeable future. However, the fore- 
ground cleaning techniques that are currently being de- 
veloped will make statistical measurement of the 21cm 
fluctuations through the power spectrum achievable with 
upcoming telescopes, and several projects are in develop- 
ment with this goal in mind, including the 21 Centimeter 
Array (21CMA; formerly PAST), the Murchison Wide- 
field Array (MWA), and the Low Frequency Array (LO- 

FAR) m. 

Complementary to power spectrum measurements will 
be the ongoing observational efforts to track the sky- 
averaged 21cm brightness temperature as a function of 
frequency (redshift). While foregrounds will make an 
absolute measurement of the brightness temperature dif- 
ficult, strong gradients over small frequency ranges will 
be detectable with the first generation of experiments 
[sgj . The hope is to use this technique to discover the 
transition between a neutral and an ionized IGM at the 
epoch of reionization, but if PBH evaporation occurs sud- 
denly enough, the sharp increase in brightness tempera- 
ture may also be a detectable signal. 

The EDGES experiment, currently under construction, 
will use a single antenna at a radio-quiet site to measure 
the sky-averaged brightness temperature and search for 
gradients of around 1 mK/MHz, with a frequency resolu- 
tion of between 0.1 and 1 MHz, in a range corresponding 
to the 21cm signal at z < 30 [H, H^]. In our example 
PBH scenarios, the brightness-temperature peak from 
sudden evaporation would be somewhat below this de- 
tectability threshold, but future experiments along these 



lines may have the power to strongly constrain sharp 
PBH mass functions through the non-detection of peaks 
in the 21cm brightness temperature. 

VII. DISCUSSION 

We have shown that future observations of the power 
spectrum of high redshift neutral hydrogen 21cm fluctu- 
ations have the potential to place strong constraints on 
the injection of energy into the intergalactic medium in 
the Dark Ages due to exotic physics, specifically in the 
case of evaporating primordial black holes. The contribu- 
tion of PBHs to the mass budget of the universe can po- 
tentially be tightly constrained in a range of masses from 
5 X 10^° kg to 10^^ kg, corresponding to PBHs whose final 
evaporation would occur during the Dark Ages when the 
universe is dominated by uncoUapsed neutral hydrogen 
gas. We estimate being able to constrain SIpbh ^ 10^^^ 
at MpBH ~ 5 X 10^° kg, and to improve upon current 
constraints derived from the 7-ray background for simi- 
lar masses. 

In order to estimate the effect of PBH evaporation on 
the IGM, we developed an extensive computer code for 
calculating the evolution and Hawking radiation of the 
PBHs, the transparency and absorption properties of the 
IGM, and the evolution of the 21cm signal both in terms 
of the sky-averaged brightness temperature and the fluc- 
tuation power spectrum. 

Although the use of 21cm observations for cosmology 
in the Dark Ages will be technically challenging, the po- 
tential payoff for a wide range of astrophysical and cos- 
mological applications easily justifies the current surge of 
interest in the field. 
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APPENDIX A: CONVENTIONS 

In this section, we describe some notational conven- 
tions and definitions used in this work, and how they 
compare to related works cited herein. 

The main source of differences among recombination 
calculations is the definition of the helium fraction, jue- 
The RECFAST code which we use for the highest 
redshifts in our calculation, defines /i^e as the helium-to- 
hydrogen number ratio: 

ijRecfast _ "^He _ Yp / A 1 ^ 

JHe = — Tn ^7\ ^ ' 

UH 4(1 -Ip) 
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where Yp is the primordial hehum abundance (mass ra- 
tio), defined by pHe/{pH + PHe), which we take to be 
0.25. Therefore, f^f^^^^ = 1/12. An alternative def- 
inition of the helium fraction is the ratio between the 
number density of helium and the total number density 
of nuclei. This definition is used by 5, hereafter FOP], 
from which we derive late-time IGM evolution equations. 
The numerical difference between the helium fraction de- 
fined this way and that of Recfast is small but potentially 
significant: 



FOP 



JHe 



riHe 



riH + nHe 



1/13. 



(A2) 



Another possible source of confusion is the definition 
of the energy injection rate. A simple option is to define 
it as the amount of energy injected per cubic centimeter, 
but it is often convenient to define it instead in terms of 
the amount of energy injected per baryon or per hydrogen 
nucleus. Investigating the energy injection due to dark 
matter annihilation in the Dark Ages, the authors of (29l . 
hereafter PF] define an energy injection rate 



where Q is the comoving energy injection rate per cubic 
centimeter, and rfi^ is the comoving number density of 
hydrogen nuclei. In a similar calculation, the authors of 
[5| implicitly define the energy injection rate due to dark 
matter annihilation or decay as 



,FOP 



Qln 



(A4) 



where is the comoving number density of baryons. 
Therefore, to compare the two calculations, one must 
keep in mind the conversion factor between the two en- 
ergy injection rates: 



,PF 



Recfast N^FOP 



(A5) 



In our calculations, we adopt the energy injection rate 
convention of [5] (equation IA4| ) and the helium fraction 
convention in RECFAST [10| (equation |Alj ). so for en- 
ergy injection rates calculated outside of Recfast and fed 
into the Recfast evolution equations, we use the conver- 

Recfast \ 



(A3) sion e ^ (1 + Affjf )e to ensure consistency 



[1] R. Barkana and A. Loeb, "The Physics and Early History 
of the Intergalactic Medium," Rept. Prog. Phys. 70, 627 
(2007) [arXiv:astro-ph/06ll54T] . 

[2] S. Furlanetto, S. P. Oh and F. Briggs, "Cosmology 
at Low Frequencies: The 21 cm Transition and the 
High-Redshift Universe," Phys. Rept. 433, 181 (2006) 
(arXiv:astro-ph/0608032 . 

[3] J. R. Pritchard and A. Loeb, arXiv:0802.2102 [astro-ph]. 

[4] A. Loeb and M. Zaidarriaga, "Measuring the small- 
scale power spectrum of cosmic density fluctuations 
through 21-cm tomography prior to the epoch of struc- 
ture formation," Phys. Rev. Lett. 92, 211301 (2004) 
arXiv:astro-ph/0312134 . 

[5] S. Furlanetto, S. P. Oh and E. Pierpaoli, "The Effects 
of Dark Matter Decay and Annihilation on the High- 
Redshift 21 cm Background," Phys. Rev. D 74, 103502 
(2006) arXiv:astro-ph/0608385 . 

[6] M. Ricotti, J. P. Ostriker and K. J. Mack, 
| arXiv:0709.0524| [astro-ph]. 

[7] S. W. Hawking, "Black hole explosions," Nature 248, 30 
(1974). 

[8] S. W. Hawking, "Particle Creation By Black Holes," 
Commun. Math. Phys. 43, 199 (1975) [Erratum-ibid. 46, 
206 (1976)]. 

[9] D. N. Page, "Particle Emission Rates From A Black 

Hole: Massless Particles From An Uncharged, Nonrotat- 

ing Hole," Phys. Rev. D 13, 198 (1976). 
[10] S. Seager, D. D. Sasselov and D. Scott, "A 

New Calculation of the Recombination Epoch," 

arXiv:astro-ph /9909275 ' 
[11] P. Sreekumar et al. [EGRET Collaboration], Astrophys. 

J. 494, 523 (1998) arXiv:astro-ph/9709257 . 
[12] B. J. Carr and J. H. MacGibbon, Phys. Rept. 307, 141 

(1998). 



[13] B. J. Carr, In the Proceedings of 22nd Texas Sympo- 
sium on Relativistic Astrophysics at Stanford Univer- 
sity, Stanford, California, 13-17 Dec 2004, PP 0204 
farXiv:astro-ph/050403£. 

[14] Ya. B. Zel'Dovich, I. D. Novikov, Soviet Astronomy, 10, 
602 (1967) 

[15] S. Hawking, Mon. Not. Roy Astron. Soc. 152, 75 (1971). 

[16] B. J. Carr, Astrophys. J. 201 (1975) 1. 

[17] D. N. Spergel et al. [WMA P Collaborat ion], Astrophys. 

J. Suppl. 170, 377 (2007) ; arXivlastra-p h/0603449 . 

[18] J. Dunkley et al. [WMAP Collaboration], 

larXiv:0 803.0586 [astro-ph]. 
[19] K. Jedamzik, Phys. Rev. D 55, 5871 (1997) 

arXiv:astro-ph/9605152 . 
[20] R. R. CaldweU and P. Casper, Phys. Rev. D 53, 3002 

(1996) arXiv:gr- qc/9509012 . 
[21] J. Garri ga and M . Sakellariadou, Phys. Rev. D 48, 2502 

(1993) arXiv:hep-th/9303024 . 
[22] S. W. Hawking, Phys. Lett. B 231, 237 (1989). 
[23] A. Polnarev and R. Zembowicz, Phys. Rev. D 43, 1106 

(1991). 

[24] J. H. MacGibbon, R. H. Brandenberger and 
U. F. Wichoski, Phys . Rev. D 57, 2158 (1998) 
[arXiv:astro-ph/9707146| . 

[25] V. A7Berezin," V. A. Kuzmin and I. I. Tkachev, Phys. 
Lett. B 120, 91 (1983). 

[26] J. Ipser and P. Sikivie, Phys. Rev. D 30, 712 (1984). 

[27] S. A. Teukolsky, "Perturbations of a rotating black hole. 
1. Fundamental equations for gravitational electromag- 
netic, and neutrino field perturbations," Astrophys. J. 
185, 635 (1973). 

[28] S. A. Teukolsky and W. H. Press, "Perturbations Of A 
Rotating Black Hole. lii - Interaction Of The Hole With 
Gravitational And Electromagnet Ic Radiation," Astro- 



15 



phys. J. 193, 443 (1974). 
[29] N. Padmanabhan and D. P. Finkbeiner, Phys. Rev. D 

72, 023508 (2005) arXiv:astro-ph/0503486 . 
[30] D. Pequignot, P. Petitjean and C. Boisson, "Total and ef- 
fective radiative recombination coefficients," Astronomy 

and Astropliysics 251 680 (1991). 
[31] Edward M. Purcell and George B. Field, "Influence of 

collisions upon population of hyperfine states," Ap. J. 

124 542 (1956). 
[32] S. A. Wouthuysen, "On the excitation mechanism of the 

21-cm (radio-frequency) interstellar hydrogen emission 

line," Astron. J. 57 31 (1956). 
[33] G. B. Field, "Excitation of the Hydrogen 21-CM Line," 

Proc. I. R. E. 46 240 (1958). 
[34] C. M. Hirata, Mon. Not. Roy. Astron. Soc. 367, 259 

(2006) arXiv:astro-ph/0507102 . 

[35] Xuelei Chen and Jordi Miralda-Escude, "The Spin- 
ndashKinetic Temperature Coupling and the Heating 
Rate due to Ly-alpha Scattering before Reionization: 
Predictions for 21 Centimeter Emission and Absorption," 
Ap. J. 602 1 (2004). 

[36] A. C. Allison and A. Dalgarno, "Spin Change in Colli- 
sions of Hydrogen Atoms," Ap. J. 158 423 (1969). 

[37] B. Zygelman, "Hyperfine Level-changing Collisions of 
Hydrogen Atoms and Tomography of the Dark Age Uni- 
verse," Ap. J 622 1356 (2005). 

[38] Piero Madau, Avery Meiksin, and Martin J. Rees, "21 
Centimeter Tomography of the Intergalactic Medium at 
High Redshift," Ap. J. 475 429 (1997). 

[39] J.R. Pritchard and S.R. Furlanetto, "21-cm fluctua- 
tions from inhomogeneous X-ray heating before reion- 
ization," Mon. Not. R. Astron. Soc. 376, 1680 (2007) 
arXiv:astro-ph/0607234 . 

[40] A. Lewis and A. Challinor, "The 21cm angular-power 
spectrum from the dark ages," Phys. Rev. D. 76, 083005 

(2007) ar Xiv:astro-ph/0 702600 . 

[41] Andrzej A. Zdziarski and Roland Svensson, "Absorption 
of X-rays and gamma rays at cosmological distances," 
Ap. J. 344 551 (1989). 

[42] X. L. Chen and M. Kamionkowski, "Particle decays dur- 
ing the cosmic dark ages," Phys. Rev. D 70, 043502 
(2004) arXiv:astro-ph/0310473 . 

[43] J. Michael ShuU and Michael E. van Steenberg, "X-ray 



secondary heating and ionisation in quasar emission-line 

clouds," Ap. J. 298 268 (1985). 
[44] K. Maeda, Classical and Quantum Gravity 3, 233 (1986). 
[45] M. J. Bowick, S. B. Giddings, J. A. Harvey, 

G. T. Horowitz and A. Strominger, Phys. Rev. Lett. 61, 

2823 (1988). 

[46] S. Coleman, J. PreskiU and F. Wilczek, Mod. Phys. Lett. 

A. 6, 1631 (1991). 
[47] M. Kavic, J. H. Simonetti, S. E. Cutchin, S. W. EUingson 

and C. D. Patterson, arXiv:0801.4023 [astro-ph]. 
[48] S. Furl anetto, Mon. Not. Roy. Astron. Soc. 371, 867 

(2006) [arXiv:astr(>-"p h/0604040 . 
[49] C. L. Carilh, New Astron. Rev. 50, 162 (2006) 

[arXiv:astro-ph/0509055 . 
[50] Y. Mao, M. Tegmark, M. McQuinn, M. Zaldarriaga and 

O. Zahn, arXiv:0802.1710 [astro-ph]. 
[51] R. T. Schilizzi, P. Alexander, J. M. Cordes, P. E. Dewd- 

ney, R. D. Ekers, A. J. Faulkner, B. M. Gaensler, 

J. L. Jonas and K. I. Kellermann, "Preliminary spec- 
ifications for the Square Kilometre Array," Square 

Kilometre Array collaboration document (2007) 
http:/ /www. skatelescope.org/pages/page_skares. htm . 
[52] B. J. Carr, J. H. Gilbert a nd J. E. Lidsey, Phys. Rev. D 

50, 4853 (1994) [arXiv:ast ro-ph/9405027 . 
[53] S. P. Oh and K. J. Mack, Mon. Not. Roy. Astron. Soc. 

346, 871 (2003) arXiv:astro-ph/0302099 . 
[54] T. Di Matteo, R. Perna, T. Abel and M. J. Rees, 

arXiv:astro-p h/010 9241 
[55] P. A. Shaver , R? AT W indhorst , P. Madau and A. G. de 

Bruyn, arXiv:astro-ph/990 1320. 
[56] X. M. Wang, M. Tegmark, M. Santos and L. Knox, As- 

trophys. J. 650, 529 (2006) arXiv: astro -ph /050108I] . 
[57] C. L. Carilh, J. N. Hewitt and A. Loeb, 

arXiv:astro-ph/0702070 
[58] J. Lazio, R. J. Macdowall, J. Burns, L. Demaio, 

D. L. Jones and K. W. Weiler, |arXiv:astro^ph /0701770| 
[59] J. D. Bowman, A. E. E. Rogers and J. N. Hewitt, 

larXiv:0710.254T] [astro-ph]. 
[60] J. D. Bowman, private communication. 
[61] The Recfast code [lol | uses the alternative convention 

/he = UHe/nH; we correct for this discrepancy in our 

calculations. 



